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ABSTRACT 

The velocity dispersion profile in globular clusters (GCs) is explained here without 
having to rely on dark matter or a modification of Newtonian dynamics (MOND). 
The flattening of the velocity dispersion at large radii in certain Milky Way GCs, 
or lack thereof, is explained by recourse to the stability of the three-body problem 
in Newtonian dynamics. The previous paper in this series determined an analytical 
formula for the transition radius between stable and unstable orbits for a star in a 
globular cluster. This stability boundary is used here to predict where the velocity 
dispersion profile is expected to flatten in GCs, given known orbital parameters of 
the GC-galaxy orbit. Published observational data for the velocity dispersion as a 
function of radius of 8 Milky Way globular clusters with approximately known orbital 
parameters has been used here. 

We find that the stability radius predicts where the velocity flattens as well as, 
and in some cases better than, the MOND models do. In particular it is noted that the 
cluster masses in the MOND models are often incompatible with the velocity dispersion 
requirements in the inner region of the cluster. However when the stability boundary is 
determined from the stability in the three-body problem using Newtonian dynamics, 
it is found to fit the observations well, and there is no need for any MOND type 
theories to explain the observations. A warning is also made against fitting the velocity 
dispersion of globular clusters without knowing their orbital parameters, particularly 
for clusters with strong evidence of ongoing tidal disruption, such as NGC 288. 
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1 INTRODUCTION 

Globular clusters are thought to contain no dark matter, 
either due to their formati on or to the evaporation of low 
mass dark matter particles (jBaumgardt fc M icskc 2008J). In 
the case of subsequent evaporation, dark matter may still 
be present in the outer regions of some clusters. However, 
for most clusters dark matter can safely be assumed to be 
absent from the cluster. This means that the cluster density 
and velocity dispersion profiles should fall off as predicted 
by a Newtonian dynamics model of a cluster in equilibrium. 

The recent discovery of flattening in the velocity dis- 
persi on profile for the ga lactic globular clusters fi Cen and 
M15 (Scar pa et al.ll2003r ). where dark matter is thought not 
to exist, has generated a lot of excitement with some stud- 
ies claiming that this is direct evidence for a breakdown 
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of Newton's laws at low accelerations. Explanations for the 
observed deviation broadly fit into three categories; tidal in- 
teractions, dark matter or a modified gravity theory (e.g. 
MOND). All of these theories can produce a flattening of 
the velocity dispersion profile for large radii, in the case of 
MON D this occurs at the critical acceleration of 1.2 x 1CP 10 
m/s 2 (|Milgromlll983l b 

Distant GCs have long bee n considered as a test 
of different gravity theories (e.g. iBaumgardt et al.l 120051 : 
iMoffat k, TothlfeoOct ). provided that the velocity dispersion 
profile is suffici ently resolved by b eing based on 30-80 stars 
in each cluster ( Haghi et al.ll201lT ). At present the statistics 
are not quite suffi cient to rule out M OND, for example in 
the case of Pal 14 jGentile et al.ll2010l '). 

A distant cluster that has received a lot of recent atten- 
tion in the literature is NGC 2419. This cluster is located 
approximately 90 kpc from the galaxy and as such, tidal 
interactions with the galaxy are expected to be negligible. 
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In addition it is a very luminous cluster and contains many 
bright giants suitable for spectroscopic analysis, greatly aid- 
ing the accuracy of the velocity dispersion profile. 

The velocity dispersion profile of NGC 2419 can only 
be fit using the MOND acceleration cut- off if there is sub- 
stant ial radial anisotropy in the cluster (ISollima fc NipotH 
The effect of radial anisotropy was found to be 
stronger in a MOND cluster profile tha n their equiva- 
lent Newtonian systems l|Nipoti et al.ll201l | ). While there is 
some radial anisotropy due to rotation, Baumgardt et all 
(2009) found that there is not enough to reconcile the 
velocity dispersion with MOND. NGC 2419 is most 
likely to be a rem nant of a much larger system as evi- 
denced dynamically jBaumgardt et al.ll2009t) and chemically 
l|Di Criscienzo et al.l 120111 ). which complicates the picture 
but still fails to reconcile the observations to any MOND 
model. A more recent and comprehensive comparison be- 
tween Newtonian an d MOND models to an increased data 
set for NGC 2419 bv llbata et"afl (|201lh found that the best 
MOND fit was far worse than the best Newtonian Michie 
model by a factor of 40000. 

For closer GCs, recent observational studies have found 
that MO ND or dark matter is not required in many GCs. For 
example lLane et alJ (|201 Oh (and references therein) found 
that of the 10 GCs studied, the velocity dispersion profiles of 
all except NGC 6121 (M4) could be well fit using a Plummer 
sphere model. They found that M4 had a mass of twice the 
literature values due to tidal heating increasing the velocity 
dispersion in the outer regions. The velocity distribution of 
this cluster was also complicated by the clear signature of 
cluster rotation found in the observations. 

The premise of this paper is that the velocity dispersion 
profile is expected to flatten out in the region where orbits 
become unstable. The method used for determining the sta- 
bility boundary is summarised he re based on the previous 
paper in this series (Kennedy 20 111 , hereafter Paper I) which 
was based on the stability o f the g enera l three-body prob - 
lem as derived in iMardlind (|2008bf) and iMardlingj ((2008a). 
This will occur in the Newtonian dynamics regime without 
need for any modifications to the theory of gravity or addi- 
tional dark matter. The determination of this radius for a 
Plummer sphere is summarised in Section [2] and is given in 
more detail in Paper I. This section also covers the velocity 
dispersion as a function of radius for Newtonian dynamics 
and the radius where this profile is predicted to break down 
according to MOND. A simulated cluster model is used in 
Section to show that the velocity dispersion profile does 
indeed flatten where stars exist in unstable orbits. 

The bulk of the paper is devoted to comparing the pre- 
dicted radius of flattening velocity dispersion to 8 individ- 
ual cluster observations from the Milky Way globular clus- 
ter system (Section [SJ. These GCs are selected on the basis 
that data is available on the velocity dispersion as a func- 
tion of radius and on the velocity of the cluster, which is 
used to determine the parameters for the cluster-galaxy or- 
bit. This section includes a discussion on any discrepancies 
between our predictions using the stability boundary and 
those by MOND for particular globular clusters. Conclu- 
sions are drawn in Section [6] in particular that modified 
gravity theories are not needed to explain the observations. 



2 STABILITY BOUNDARY IN A PLUMMER 
SPHERE 

The previous paper in this series outlined a simple method 
for determining the transition between stable interior and 
unstable exterior orbits in a GC. An unstable system is 
taken to mean that a star is on an orbit that makes the total 
system unstable to the escape of one of the bodies. In the 
context of GCs this is equivalent to the eve ntual escape of a 
star from the potential well of the cluster. iKenne 13 (f201lh 
used this stability boundary to estimate the tidal radius for a 
cluster given its galactic orbital parameters of perigalacticon 
(R p ) and eccentricity (e). This radius was derived assuming 
point masses for the distant galaxy and stars orbiting within 
the cluster, which is described by a Plummer potential cen- 
tred on the cluster core. In this section we give an overview 
of the results from this paper where they are relevant to 
predicting the flattening of the velocity dispersion in GCs. 

A globular cluster is modelled usin g a Plummer sphere 
with gravitational potential given by (Bi nnev fc Tremaind 
1 19871) 



$ = 



-GM C 



r i/2 + r2 



(1) 



where G is the gravitational constant, Mc is the mass of 
the cluster, r is the radial distance, and r 1 / 2 is the observ- 
able half mass radius of the cluster for a projected velocity 
disp ersion. 

iDeionghel l| 19871) gives the velocity dispersion as a func- 
tion of cluster radius, r, for a Plummer sphere to be 
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and the central velocity dispersion (ao) is related to the 
cluster mass Mc by 
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where G has its usual value. 

Throughout this paper there are two radii of key im- 
portance for each GC. These are the radius where the ac- 
celeration acting on a star due to the cluster potential goes 
beneath the MOND limit of a = 1.2 x 10" 10 m/s 2 (r a0 ), 
and the radius of the transition between stable and unsta- 
ble orbits (r c haos, see below). Firstly the MOND radius is 
found by iterating over the acceleration ao and solving for 
r a o, which for a Plummer sphere is 



GM c r a0 
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given the known values of ri/2 and Mc- 

For the second radius we adopt a functional form for 
the tidal radius that is separable into a mass component 
and orbital eccentricity of the cluster-galaxy orbit 



n 
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where Mc and Mc are the masses of the cluster and galaxy 
respectively, e is the eccentricity of the cluster orbit around 
the galaxy, and R p is the distance of closest approach to the 
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Table 1. Coefficients for the polynomial fit to r c i laos and the 
minimum and maximum extents of the marginally chaotic zone. 
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For radii r < r m in all orbits are expected to be stable, whereas 
for r > r max they are expected to be chaotic. The resultant 
functions of the eccentricity of the cluster-galaxy orbit are shown 
in Figure [T] For de tails on how these coefficients are determined 
see lKennedvl ll201ll1 . 

galaxy, referred to as perigalacticon. This form of the tidal 
radius also allows direct comparison between the stability 
boundary and other tidal radii estimates from the literature. 

The most commonly used tidal radius estimate is for a 
star on a radial orbit and using point mass potentials for the 
cluster and galaxy. The eccentricity dependence, /(e), of the 
tidal radius is given by the classical iKind l|l962f ) result 

/(e) = fc(3 + e)- 1/3 (6) 

where the constant k ~ 0.7 was introduced by iKeenanl 
(1981) to better fit observations of the galactic globular clus- 
ters. The King radius (rt) will be used to denote the maxi- 
mum theoretical tidal radius of a GC by using Equations ((5} 
and (|6]) with k = 1 and e = 0. 

The stability boundary was determined in Paper I us- 
ing the eccentricity distribution for stars in an isolated clus- 
ter as modelled by a Plummer sphere of potential given by 
Equation [T] The transition between stable interior orbits 
and unstable exterior orbits was not a single value; rather it 
was a range of radii that depended on the eccentricity of the 
cluster-galaxy orbit. To describe this transition three values 
were used, the minimum and maximum values and an in- 
dicative value denoted by r c h aos . This value was defined as 
the lowest radius where the fraction of unstable orbits drops 
beneath 10%. 

Using the same form for the tidal radius as Equation ((5} 
with the same mass dependence then the eccentricity depen- 
dence for all three values can be fit using a Taylor series in 
e of the form 



/(e) = exp 
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where the coefficients are given in Table [T] for the indicative 
stability boundary radius (r c haos) and the minimum and 
maximum extents of the partially stable region. For radii 
r < r m i n all orbits are expected to be stable, whereas for 
r > rmax they are expected to be chaotic. The eccentricity 
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eccentricity of the cluster-galaxy orbit are shown in Figure [T] 
The functional dependence on the eccentricity was de- 
rived using a distribution of stellar orbital eccentricities and 
was effectively independent of the cluster mass. The assump- 
tion of a point mass potential for the galaxy for distant clus- 




Figure 1. The eccentricity dependence for the transition from 
stable inner orbits to unstable exterior orbits. The form of the 
eccentricity dependence is given by Equation [7] with coefficients 
given in Table [T] and the corresponding radius is given by Equa- 
tionO The indicative value for the stability boundary (associated 
with r c haos) is shown as black, while the minimum and maximum 
extents of the partially stable region are the red and green curves 
respectively. 



ter orbits is justified for distant GCs since the mass inside 
6.4 kpc is consistent with that of a point mass of roughly 
10"M fl , based o n the orbits of NGC 2419 and NGC 7006 
(Bcllazzim 120041 ). However all of the GC orbits examined 
here have perigalacticon distances (R p ) inside 6 kpc. The 
procedure adopted in Section [4] determines an effective Ke- 
plerian orbit by setting R p = R p (min) and getting the ec- 
centricity from the apogalacticon R a — Rajnax from mul- 
tiple passages. This means that the orbital parameters are 
themselves an approximation of a Keplerian orbit to a more 
complicated orbit. In this context the error in the GC-galaxy 
orbit due to modelling the galaxy as a point mass is negligi- 
ble. The overall effect of the uncertainties in the GC-galaxy 

orbit is discussed in Se ction [5] 

As pointed out bv lFukushige fc Heggid (|2000h the tidal 
radius does not act as an instant remover of stars. They 
found for GCs on circular orbits that the escape timescales 
for stars beyond the tidal radii could be long enough to al- 
low some stars to stay in this region indefinitely. This means 
that the predicted tidal radius will be a lower limit as stars 
outside this can still be close to the GC while remaining 
formally unbound. Stars on such unstable orbits will be ob- 
servable inside the King radius, but are expected to have 
very different velocities compared to an equilibrium cluster. 
The effect of stars on these orbits on the observable velocity 
dispersion profile is discussed in the next section. 



3 VELOCITY DISPERSION PROFILE IN A 
SIMULATED CLUSTER 

It is expected that stars on chaotic orbits at large distances 
from the cluster centre lead to a different velocity profile 
than an equilibrium model (such as Equation [2}. This sec- 
tion uses a simple cluster model to estimate what effect stars 
on chaotic orbits in the outer regions of a globular cluster 
has on the velocity dispersion profile. 

A globular cluster is modelled using a Plummer sphere 
with gravitational potential given by Equation Q centred 



4 Gareth F. Kennedy 



on a single particle that in turn orbits the galaxy. To examine 
the effect of the tidal field, without prescribing the force 
as a function of time, a free moving point mass particle of 
mass 10 11 Mq is used. The stars in the outer regions of the 
cluster are modelled as point masses of IMq that respond 
to the movement of the cluster centre particle of 1O 6 M0, 
but not to each other. This same model was used in Paper I, 
without the inclusion of the galactic tidal field, to determine 
the distribution of orbital eccentricities for stars inside the 
cluster in finding the stability boundary. 

By removing the mutual interactions between particles 
the cluster model does not include two-body relaxation. This 
is not a problem for the outer regions of interest here since 
typical relaxation timescales are greater than 1 Gyr at the 
half mass radius and scale as approximately r 1//2 . Therefore 
the relaxation timescale in the outer cluster regions is close 
to the age of the cluster itself and the region is collisionless. 
There is a form of relaxation due to the cluster centre par- 
ticle moving in response to the distribution of cluster halo 
stars. The wandering of the cluster centre moves the Plum- 
mer potential which accelerates the loss of stars over the 
tidal radius. From extensive simulations it was found that 
the wandering core is not a problem until the cluster has lost 
most of its original mass. At this point the cluster is strongly 
disrupted and its velocity dispersion cannot be adequately 
described by a simple equilibrium model. 

The primary interest here is to observe deviations from 
the velocity dispersion for a cluster in equilibrium, given 
by Equation ([2]), thus the distribution of stars is artificially 
biased to provide better resolution in the outer regions. A 
uniform distribution in radius between ri/2 and the King 
radius rt is chosen and sampled with 10 3 stars. The region 
inside ri/j is modelled by an unchanging potential rather 
than particles since the velocity distribution inside the half 
mass radius remai ns isotropic and theref ore does not need 
to be resolved (e.g. iGiersz fc Heg gie 1994|). All particles arc 
numerically integrated for a total of 30 cluster-galaxy or- 
bits (« 4.4 Gyr) using a collisionless N-bo dy code with a 
Bulirsch-Stoer integrator (|Press et al.lll98r3 ) that conserves 
energy and angular m omentum. Furth er details about the 
model are described in I Kenned vl |2008t ). 

Many models were run to sample across a range of 
cluster-galaxy orbital parameters; here we focus on one 
model where the stability boundary was well resolved. The 
system parameters are R p = lOOOr - !^, e out = 0.3 and 
ri/2 = Ipc, M c = 10 6 M© and M G = 1O U M . The velocity 
dispersion profile against radial distance to the cluster cen- 
tre after 10 (black), 20 (red) and 30 (green) cluster-galaxy 
orbital periods are shown in Figure [2] All time snapshots 
are taken at the apogalacticon to reduce any spurious per- 
turbations to the halo stars from perigalacticon passages. 
The transition between stable and unstable orbits (r c haos) 
is shown as a solid vertical line, with the range (r min / max ) 
shaded to show the size of the transition. The radius for 
which the acceleration is equal to the MOND value, given 
by Equation [4l is shown as a dashed vertical line. 

As is evident from Figure [5] the stability boundary is a 
good predictor for the flattening of the velocity dispersion. 
Coincidently the MOND radius for this cluster also occurs 
before the flattening despite the purely Newtonian dynamics 
included in the simulation. The occurrence of flattening in 
the velocity dispersion and its accurate prediction by the 
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Figure 2. The velocity dispersion profile against radius from a 
collisionless N-body cluster simulation composed of 1000 parti- 
cles. The velocity dispersion profile after 10 (black), 20 (red) and 
30 (green) cluster-galaxy orbits is shown. The transition from 
stable inner to unstable outer orbits is shaded and the indicative 
radius (r c haos) is shown as a vertical line. The radius for which 
the acceleration is equal to the MOND value is shown as a dashed 
vertical line. Note that the velocity dispersion begins to flatten in 
the outer unstable region. 



stability boundary provides confidence applying this method 
to observations of GCs in the Milky Way system in Section[5] 



4 DETERMINATION OF GC-GALAXY 
ORBITS 

Before comparing the velocity dispersion profiles of real clus- 
ters to the radii estimates in Section[2]the orbital parameters 
of the cluster-galaxy orbit are required. In this section the 
galactic orbits of 8 GCs are determined based on the ob- 
served velocities for each cluster and integrating backwards 
in time through a realistic galactic potential. The perigalac- 
ticon (R p ) and eccentricity (e) can then be obtained from 
the integrated orbits for each GC. 

Recent velocity and distance values for each cluster 
(references in Table [2]) are used in con j unctio n with the 
galactic potential used bv lFellhauer et all (120071 ). This grav- 
itational potential consist s of a Miyamoto-Nagai potential 
l|Mivamoto fc Nagailll975h combined with a logarithmic po- 
tential. The total galactic potential $ is given as a sum of 
the galactic halo $h, disc $d and bulge <I>6 potentials by 
l|Fellnauer et al.ll2007l ) 



$(2, y, z) = $ h (r) + $ d (R, z) + $,,(r) 



where 

r — \J x 2 + y 2 + z 2 



R = yjx 2 + y 2 



(8) 



(9) 



and x, y and z are galactic coordinates with units in kpc. In 
this coordinate system the Sun is located at (-8,0,0) at which 
a particle with velocity directed in the positive y direction is 
moving in the direction of Galactic rotation. It follows that 
z points in the direction of the northern galactic pole. The 
gravitational potentials for the galactic halo <&h, disc $d and 
bulge $(, are 
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where a = 12.0 kpc, b — 6.5 kpc, c = 0.26 kpc, d — 0.7 
kpc, v — 181 km/s and G is the gravitational constant. 
The masses of the galactic disc and bulge are Md = 10 11 
M Q and M b = 3.5 x 10 10 M Q respectively. The equations of 
motion for a cluster moving in this potential are given by 



(13) 



For globular clusters with velocity data given in the liter- 
ature, Equation (|13[) is integrated back through time for 
approximately 10 cluster-galaxy orbits, using the Bulirsch - 
Stoer numerical integration method |Press et al.l [l986). 
Then r m in and r max are measured and used to approximate 
the effective Keplerian orbital eccentricity and perigalacti- 
con. 

The initial conditions and resulting orbital parameters 
from these cluster-galaxy orbit simulations are summarised 
in Table [2] This table lists the observed cluster mass (Mo), 
the distance to the galactic centre (Rg), orbital parame- 
ters (Rp, e), half- mass radius (ri/ 2 ), the region where the 
velocity dispersion flattens (r/iat), the theoretical estimates 
where this will occur (r chaos and r a o), the tidal radius given 
by Equation © and finally the references for orbital and 
velocity dispersion observational data. The orbital param- 
eters R p , e and P for the G Cs in Table [ 2] gene rally agree 



with the parameters found bv lAllen et al. (2006) who used 
a different galactic potential. 

Note that there are two sets of masses (Mc) for each 
cluster in Tabled The first of t hese sets contains th e MOND 
paper cluster masses quoted in lScarpa et alj l|201lf ) (and ref- 
erences therein) and fitted such that r a o is close to rfi a t- The 
second mass set is determined using the central velocity dis- 
persion for each cluster from a fit to all data points within 
r a o, where the cluster is in the Newtonian regime and Equa- 
tions ([2| and ([3]) are valid. 

Most of the uncertainty in the orbital parameters (R p 
and e) is due to the observations of the cluster velocity com- 
ponents directed tangential to our line of sight having large 
uncertainties. Such measurements are difficult to obtain and 
rely on comparing the current positions of globular clusters 
on the sky with older position data found on photographic 
plates. The motion of each cluster is then determined from 
the change in positions (in milli-arcseconds) over time, with 
the positions c alibrated using "fixe d" distant objects, such 
as quasars (e.g. iDinescu et all 19971 ). The error is also ampli- 
fied by uncertainties in the distance to the cluster, which also 
affects the half- mass radius measurement. From the point of 
view of the orbital parameters in Table [2] these uncertain- 
ties are most likely to translate to an error in eccentric- 
ity if Rp/Rc ~ 1 and in eccentricity and perigalacticon if 
Rp « Rg. It is proposed here that the error in the cluster- 
galaxy orbital parameters is the reason for discrepancies in 
the velocity dispersion profile for particular clusters in the 
next section. 



5 GC VELOCITY DISPERSION PROFILES 

This study is limited to clusters for which online data is 
available for the velocity dispersion as a function of ra- 
dius where flattening has been reported. We are further 
limited to clusters where the orbital parameters can be de- 
termined (previous section). The galactic globular clusters 
fitting these requirements are listed in Table [2] 

Observational data for the velocity dispersion pro- 
files for each cluster are given in Table [2] along with the 
data sources for each clust er. Most veloc i ty dis persion data 
available online was from IScarpa et al.l I^OHl - ) (and refer- 
ences therein) who chose cluster masses consistent with the 
MOND radius fitting the flattened region. These masses are 
collectively referred to as the MOND cluster masses and 
are listed first in Table [2] for each cluster. Distances to the 
galact ic centre (Rg) for all clusters were taken from lHarrisI 
l| 19961 ). References for the orbital parameters are listed in 
the table as numbers rather than letters for each cluster. 

Two special cases for the observed velocity dispersion 
data exist in this sample. Firstly, the velocity dispersion 
data for NGC 707 8 are taken from the non-rotating cluster 
model favoured bv lDrukier et al.1 (1 19981 ). rather than the ro- 
tating cluster model. Secondly, th e velocity dispersion data 
for NGC 6341 l|Drukier et alj|2007| ) was not binned in radius 
in the same way as the other clusters in this sample. For this 
cluster no data was given on the minimum or maximum dis- 
tances for each bin, nor for the number of stars used in each 
bin. This means that no radial error bars are available for 
the data points for NGC 6341. 

Figures [3] and [4] show the velocity dispersion profiles for 
all clusters listed in Table[2]using the MOND cluster masses. 
The MOND radius r a o as calculated for each cluster from 
Equation ((3]) is shown in each plot as a dashed red line. For 
the MOND masses this is close to where the velocity disper- 
sion flattens out by the choice of the cluster mass Mc- The 
transition from stable to unstable orbits is shaded in grey 
and the indicative radius r c haos is shown as a solid line inside 
this region. The theoretical maximum extent of the cluster, 
the tidal radius using Equation , is shown as a dotted ver- 
tical line when it is inside the plotted range. Finally the red 
dotted curve indicates the theoretical velocity dispersion for 
the quoted cluster mass using the radial dependence given 
by Equations (|2j) and (|3} . 

The second set of masses in Table [2] uses cluster mass 
values determined by fitting the observational data points 
within r„o. The fits are calculated by Equations ((2| and ([3} 
and the error is minimised using a \ 2 test for all data with 
r < r a o- However since r a o depends on Mq which in turn 
depends on the data points used to get the fit from which 
(To is taken then Mc must be solved in an iterative manner. 
The resulting masses are listed in the table and are also 
shown next to the names of each cluster in Figures [5] and [6] 
These figures present the same data as previously shown in 
Figures [3] and 2] differing only in the cluster mass used to 
determine the stability transition region, r c h a0 s, Vt, r a o and 
the theoretical velocity dispersion profile. 

The results from Figure [3l6l can be summarised in terms 
of how close the MOND acceleration radius (r a o) or the 
transition to unstable orbits (r c haos) fits the flattening of 
the velocity dispersion. Differences between the MOND and 
chaotic orbit predictions can be described in three general 



6 Gareth F. Kennedy 

Table 2. Globular cluster parameters for 8 Milky Way globular clusters. 
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Cluster masses from fits using MOND and from the inner velocity dispersion are given in column 2. Orbital and structural parameters 
for each cluster are given in columns 3-6. The estimate for the observed flattening is given in column 7 and predict ions for th is value 

using the stability boundary and the MOND radius where ao = 1.2 X 10~ 8 cm s — 2 assuming M/L = 1 llScarpa et alj|201ll . and 
references therein) are given in columns 8 and 9 respectively. The tidal radii as calculated by Equation J5} is given in column 10. The 
final column contains the references for the orbital parameters (often the raw cluster velocitie s) and the velocity d ispersi on profiles. Th e 
refere nces a re den oted by numbers in the fo rmer case and l etters in th e later and are; A = Drukier et al. I dl998f ) . B = IScarpa et al] 

d2003h. C =IScarpa et"aD ll2004fl. D = IDrukier et al] ll2007ll. E = IScarpa et all l|2007ah . F = IScarpa et al] d2007bh. G = 

IScarpa fc Falomol d2010h. H =IScarpa et all [1201111 , 1 = lMeziane fc Colinl lll99rj) and references therein. 2 = |Pinescu et al] (Il999ll and 
references therein, 3 = |pinescu et aU^2*001*l v See in-text for more details including how the orbital parameters R p and e are determined. 



categories; estimates are consistent in cluster masses, incon- 
sistent masses, and problem clusters for both models. 

5.1 Consistent masses 

This category is defined as the set of clusters where the 
MOND mass and the fitted mass are in good agreement and 
where both r a o and r c h aos are consistent with the flattening 
of the velocity dispersion. This category includes; NGC 288, 
NGC 1904, NGC 7078 and NGC 7099. 

The velocity dispersion of NGC 288 is consistent with 
being completely flat with an average value of 2.3 ± 0.15 
km/s (|Scarpa et al.|[2007al ). This cluster has a low concen- 
tration cluster with a mass of ~ 5 x 1O 4 M0 and a half-mass 
radius of riyj = 5.76 pc. As such the acceleration due to 
the cluster potential is below the MOND limit for all values 
of radius, i.e. r a o = pc. The stability boundary and the 
tidal radius are also very close to the centre of the cluster, 
as seen in panel (a) of Figures [3] and [5] NGC 288 has a 
galactic orbit with a low perigalacticon (1.7 kpc) and high 
eccentricity (0.74) which means that its tidal radius is in- 
side the observable cluster limits; hence it is currently being 
disrupted. Additiona l velocity di s persio n observational data 
to 26 pc obtained bv lLane et al. (2010) found that a Plum- 
mer model provided a good fit to a low concentration cluster 
without significant deviation from Newtonian dynamics. 

In the case of NGC 7078 r a o is closer to the cluster 
centre than r chaos, but both radii are consistent with the 
flattening of the velocity dispersion. The observational data 
presented in the figures is b ased on the no rotation model 
favoured by IDrukier et al] |l998). They also put forward 
the explanation of the flattening in the velocity dispersion 



as being due to tidal hea ting, as expected in the models of 
I Allen fc Richstonel (|l988h . 

For the remaining clusters with consistent masses, 
NGC 1904 and NGC 7099, we find that r a0 « r chao3 . 
In other words either MOND or the occurrence of chaotic 
orbits in the outer regions can be used to explain the flatten- 
ing of the velocity dispersion. This is the general conclusion 
for all four clusters in this category. 

5.2 Inconsistent masses 

This category is defined as those where the MOND mass and 
the fitted mass are in substantial disagreement. In general 
the r a Q is a good fit to the velocity dispersion flattening 
only when the MOND masses are used, however these are 
inconsistent with the velocity dispersion profiles in the inner 
regions where it should by purely Newtonian. This category 
includes NGC 1851 and NGC 6171. 

MOND has particular difficulty fitting the flattening 
of NGC 1851 with the fitted cluster mass value of 3.6 X 
10 5 Mq, however the same is true for the stability transition 
for this cluster. For the MOND cluster mass of 1.8 X 1O 5 M 
both r a o and r c haos are consistent with the flattening of 
the velocity dispersion at about 18 pc. It is proposed that 
there is most likely an error in the orbit (R p = 5.7kpc and 
e = 0.69) for this cluster. For example changing the peri- 
galacticon of the orbit to R p = 5 kpc or increasing the ec- 
centricity both decrease r c h a0 s- However changing the orbit 
will not change the radius where the acceleration drops be- 
neath the critical MOND value. 

The MOND cluster mass given by IScarpa et aU (|2004T ) 
for NGC 6171 seems peculiarly lower than the fitted mass, 
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(b) NGC 1851; M c = 1.8 x 1O 5 M 
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Figure 3. The velocity dispersion profile for all globular clusters in this study using the cluster mass estimates given in MOND papers. 
Masses arc shown beside the cluster name and the other cluster parameters are given in Table [2] The transition from stable to unstable 
orbits is shaded in grey and the indicative radius r c f laos is shown as the solid vertical lines. The vertical red dashed line shows where the 
acceleration is equal to the MOND value. The red dotted curve indicates the theoretical velocity dispersion for the quoted cluster mass. 



particularly since the MOND radius gives a good fit to the 
flattening of the velocity dispersion in both cases. Interest- 
ingly for both masses the stability transition is interior to 
the MOND radius and as such provides a better fit to the 
flattening. 



transition to chaotic orbits than by the radius where the 
acceleration equals the MOND value. 



5.3 Problematic clusters 



The flattening of the velocity dispersion for both of This category includes the clusters which present consis- 
these clusters is found to be a better fit by the stability tency problems for both the MOND model and the chaotic 
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(a) NGC 7078; M c = 3.9 x 1O 5 M 



(b) NGC 7099; M c = 7.8 X 10 4 M Q 
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Figure 4. Figure [3] continued; MOND cluster masses used. 
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orbits model. The two clusters in this category are NGC 
5139 (w Cen) and NGC 6341 and a longer discussion of the 
results for these clusters is warranted. 

As is typical for this cluster, the velocity dispersion pro- 
file of NGC 5139 (uj Cen) presents difficulties for both 
models. In the case of MOND the radius r a o can be made 
to roughly match the observed flattening near 30 pc only by 
using a mass of Mc = 1.1 X 10 6 Mq. However this mass 
is not compatible with the velocity dispersion inside the 
cluster in the Newtonian regime, which requires a mass of 
Mc = 2.5 x 10 6 Mq. Using chaotic orbits for this cluster 
is also difficult in that the stability boundary is found to 
be well inside the cluster for the orbital parameters given 
in Tabled Two possible resolutions to this incompatibility 
might come from th e fact that NGC 5139 is likely to be 
rapidly rotating fe.g. lScarpa fc Falomoll2~010r i and/or that it 
is suffering ongoing disruption (see below). 

The final cluster examined in this study is NGC 6341 
which has a cluster mass of approximately Mc = 1.5 x 
10 s Mq by both the MOND literature and the fitting proce- 
dure used here. The MOND radius of r a o = 13 pc is signifi- 
cantly exterior to the observed velocity dispersion, flattening 
at about 8 pc, which is itself exterior to the predicted sta- 
bility boundary of r c ; laos = 4 pc. A key observation for this 
cluster is that the flattening occurs very close to the tidal 
radius itself. This implies that the cluster is being severely 
tidally disrupted which explains the flattening without need 
for MOND. 

Both of these clusters, especially NGC 6341, have evi- 
dence of ongoing tidal disruption. Note that neither the tidal 
radius nor the stability b oundary act as an instant r emover 
of stars, as pointed out bv lFukushige fc Heggiei (|2000h . They 
found for GCs on circular orbits that the escape timescales 
for stars beyond the tidal radii could be long enough to allow 
some stars to stay in this region indefinitely. For eccentric 
cluster-galaxy orbit this timescale would be shorter while 
for internal orbits within the cluster the timescale would be 
substantially longer. It seems likely that this effect is seen 
in both NGC 5139 and NGC 6341. 

An alternative explanation for these clusters is that the 
orbital parameters are in error. As mentioned previously the 
eccentricity (e) and to a lesser extent the perigalacticon (R p ) 
have large uncertainties. Changing either of these quantities 



has a dramatic effect on the estimate for r c h aoa but no effect 
on the position of r a o when only the cluster potential is con- 
sidered. For example an increase in eccentricity of Ae ~ 0.1 
will produce approximately a 10% decrease in r c h aos , for a 
fixed value of R p . While this may be sufficient to explain 
NGC 6341, it is not sufficient to explain NGC 5139 which 
may require a higher initial mass and ongoing tidal disrup- 
tion to explain the velocity dispersion. 



6 CONCLUSIONS 

It is found that MOND is not required to explain the flat- 
tening of the velocity dispersion profile in globular clusters. 
We examined eight GCs from the Milky Way GC system 
with observational data for both the velocity dispersions as 
a function of radius from the cluster centre and with data 
for the orbital parameters for the cluster-galaxy orbit. 

From all of these clusters we found that the flattening 
could be equally, if not better, predicted using the tran- 
sition from stable interior to unstable exterior orbits than 
using the radius where the acceleration equals the MOND 
value. Of the eight clusters, a minor change in the perigalac- 
ticon and eccentricity of NGC 1851, and more substantial 
changes to those for NGC 6341, would result in a much bet- 
ter fit using the stability results in the Newtonian regime 
than the MOND predictions. Since the orbital parameters 
are determined from observations of the proper motions of 
the clusters, these quantities are known to be very uncertain. 
In the case of NGC 6341 the errors in the orbital parameters 
need not be so large if the cluster is undergoing substantial 
tidal disruption, as seems likely. 

One missing feature in the analysis of the velocity dis- 
persion profiles of GCs presented here is the non-rigorous 
treatm ent of cluster rotation. For example ISzekelv et al.l 
|2007l ) studied five clusters (NGC 104, NGC 288, NGC 6218, 
NGC 7099 and NGC 6809) and found evidence for rotation 
in four of them. They also emphasise the point that to disen- 
tangle these effects from a possible breakdown of Newtonian 
dynamics would require observations of a much larger range 
of globular clusters. We would add to this that the orbital 
parameters must be known in order to do any meaningful 
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(b) NGC 1851; M c = 3.6 x 10 s M 
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Figure 5. Velocity dispersion for all globular clusters in this study using mass estimates from the fit to the <r(r) values inside of r a Q. 
Masses are shown beside the cluster name and the other cluster parameters are given in Table [2] The transition from stable to unstable 
orbits is shaded in grey and the indicative radius r c (, aos is shown as the solid vertical lines. The vertical dashed red line shows where the 
acceleration is equal to the MOND value. The red dotted curve indicates the theoretical velocity dispersion for the quoted cluster mass. 



analysis involving tidal effects on globular clusters, in par- 
ticular on the velocity dispersion. 



We find that an error in the cluster-galaxy orbital pa- 
rameters, known to be very uncertain, is more likely than 
requiring new physics to explain the velocity dispersion in 
some clusters. 
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